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Abstract 

We show that a large subclass of variograms is closed under Schur products and that some desirable stability properties, like 
the Schur product of ad hoc compositions, can be obtained under the proposed setting. We introduce new classes of kernels 
of Schoenberg-Levy type and show some important properties of eventually constant, radially symmetric variograms. In 
particular, we characterize eventually constant variograms in terms of their permissibility in Euclidean spaces of arbitrary 
high dimension. 
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1 Introduction 

Positive and conditionally negative definite functions on groups or semigroups have a long history and they appear 
in many applications in probability theory, operator theory, potential theory, moment problems and various other 
areas. They enter as an important chapter in all treatments of harmonic analysis, and their origins can be traced 
back to papers by Caratheodory, Herglotz, Bernstein and Matthias, culminating in Bochner's theorem from 1932, 
see the survey in Sasvari [22]. Schoenberg's theorem explains the possibility of constructing rotationally invariant 
positive definite and conditionally negative definite functions on Euclidean spaces via completely monotone func- 
tions and Bernstein functions. Positive and conditionally negative definite functions are a cornerstone of spatial 
statistics where they are, respectively, known as covariances (or kernels) and variograms. The theory of random 
fields, started in the in the 40s with the early works of Kolmogorov (see [6] and references quoted there) and 
further developed by Gandin [8] and Matheron [18j among others, is based on the specification of these classes. 
In particular, the kriging predictor, viz. the best linear unbiased predictor, depends exclusively on the underlying 
covariance or variogram and we refer to the tour de force in Stein [26j for a rigorous assessment of this framework. 

Let {Z{^), ^ G K'^} be a weakly stationary and Gaussian random field (RF for short). The associated covariance 
function C : R'' M is positive definite, that is for any finite collection of points G K'* the matrix 

(C(^i — is positive definite, i.e. 



Thus, a function C : K'' R is positive definite if and only if there exists a weakly stationary and Gaussian RF 
having C(-) as covariance function. If C(-) is radially symmetric, the associated Gaussian RF is called isotropic. 

It is well known that the family of covariance functions is a convex cone which is closed under addition, products, 
pointwise convergence and scale mixtures; for these basic facts the reader is referred to standard textbooks of 
geostatistics such as Chiles and Delfiner [6]. 

A variogram 7 : R'' ^ C is the variance of the increments of an intrinsically stationary RF, that is, for any two 
points £,1,^2 G R'^, Var(Z(^i) — -^(^2)) ■— 7(^1 — ^2)- Variograms are conditionally negative definite, that is for any 
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for all ai , 02 , . . . , a„ S C : 




1 



finite collection of points {■^i}"^! G K'^ we have 

n n 

for all ai, . . . , a„ e C such that a, = : ai7(Ci ~ — 0- (1) 

i—l *,J — 1 

There is a close relationship between variograms 7 and stationary covariance functions C. The elementary estimate 
|C(C)I — ^(0) =: VarZ shows that stationary covariance functions are necessarily bounded; in particular, 7(^) := 
C(0) — C(^) is a variogram. Indeed, variograms may be unbounded, as in the case of fractional Brownian motion. 
If, however, the variogram is bounded, then it is necessarily of the form C(0) — C(^), ^ G R'' for some stationary 
covariance function C(-). 

The terminology concerning positive and conditional negative definiteness is not uniform throughout the lit- 
erature; it depends very much on the mathematical context or the application in science. Christakos [7] and 
many other applied scientists use the notion of permissibility for both concepts: positive and conditional negative 
definiteness. We will use both conventions alongside whenever no confusion can arise. 

If C is compactly supported, i.e. if it is zero outside some compact set K C K'', the associated variogram 
7(^) :~ C'(O) — C(^) will be called from now on eventually constant. Compactly supported covariance functions 
have been well understood since the early 90's with the works of Gaspari and Cohn [9] and Gneiting [12] and they 
are particularly important in recent literature oriented to the analysis of massive spatial datasets. Thus, eventually 
constant variograms represent an important tool whose properties have not been assessed yet, and these facts 
constitute an important motivation for the present work. 

In this paper we are mainly interested in radially symmetric covariances and variograms. This means that the 
associated Gaussian RF is weakly or intrinsically stationary and isotropic. Isotropy and stationarity are independent 
assumptions but we will assume both to keep things simple. An isotropic covariance function, rescaled by its value 
at the origin, is the characteristic function of a radially symmetric random vector on the sphere of R'^. This class 
of covariances is well understood and we refer to Gneiting \W\ \TT\ and the references given there for an extensive 
survey of this topic. Much less is known about variograms. For instance, it is common knowledge that the class of 
variograms is a convex cone which is closed in the weak topology of pointwise convergence, but the Schur product 
of two variograms is not necessarily a variogram. This is a point that deserves a thorough discussion in the light of 
a recent beautiful result in Ma [iTt Theorem 3(i)] where a simple permissibility condition is given for the product 
of two exponential variograms composed with a homogeneous function. 

We shall also give a general answer to this question as well as a complete characterization of those variograms 
whose Schur product preserves conditional negative definiteness. We shall then focus on other challenging problems 
related to ad hoc compositions of variograms as well as to quasi arithmetic compositions of them. 

The use of kernels of Schoenberg-Levy type has been persistently emphasized by both old and recent literature. 



2 



In this paper we give new forms of kernels of this type that may be appeaUng for modelling in spatial statistics. 

Another crucial problem faced in this paper regards the potential trade-off between, on the one hand, the 
computational advantages induced by the use of compactly supported kernels and, on the other hand, the fact that 
compactly supported kernels can be positive definite only on finite dimensional spaces, by a striking and beautiful 
result due to Wendland [29]. We consider this problem from the point of view of variograms; this makes sense 
since variograms, which are possibly unbounded, represent a larger class than covariance functions. As mentioned 
before, compactly supported covariances correspond to eventually constant variograms. We will characterize all 
eventually constant radially symmetric variograms in dimensions d > 2 and give some characterizations of these 
classes in terms of their eventual periodicity. 

The paper is organized as follows: Section 2 contains the basic material required for a self-contained exposition 
and for understanding the technical proofs of our statements. Section 3 assesses new stability properties of the 
variogram class, while Section 4 is dedicated to kernels of Schoenberg-Levy type. Section 5 discusses the relation 
between eventually constant variograms and their permissibility in arbitrarily high dimensional spaces. We also 
assess some other important properties of this class. 

2 Pick-Nevanlinna functions and complete monotonicity 

This section is mainly expository and we collect here some basic material needed later on. We will frequently use 
the following characterization of variograms, for which a proof can be found in Berg and Forst [S] Proposition 7.5]. 

Theorem 1 . A function j : M."^ M. is a variogram if, and only if, the following three conditions are satisfied 

i) 7(0) > 0; 

ii) 7(0 = 7(-0; 

iii) 7 is conditionally negative definite, i.e. equation ^ holds for a/Z ^i, . . . , ^„ € R''. 

Often Polya's Theorem, cf. [31 Theorem 5.4] comes in handy if one wants to construct concrete examples of 
variograms. 

Theorem 2. A continuous function : M ^ [0,oo) which is even, i.e. 4>{x) = a;), decreasing and convex on 
the interval (0,oo), is positive definite. 

Clearly, 0(0) — (j){x) is increasing, concave and a variogram, e.g. [51 Corollary 7.7]. 

Recall that a function / :]0, cx)[-^ R is called completely monotone, if it is arbitrarily often differentiable and 

(_l)«/("'(a;) > for a; > 0, n = 0, 1, . . . . 
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By Bernstein's theorem the set CM of completely monotone functions coincides with the set of Laplace transforms 
of positive measures /i on [0,oo[, i.e. 



where we only require that c""^* is /i-integrable for any x > 0. CM is a convex cone which is closed under addition, 
multiplication and pointwise convergence. 

Definition 3. A function / : ]0, oo[ ^ M is called a Stieltjes function, if it is of the form 

m-a+f 'Jf- (2) 
where a > 0, and ^ is a positive measure on [0, cx)[ such that j^^ + t)^^ dfi(t) < oo. 

The set § of Stieltjes functions is a convex cone such that S C CM. Since Stieltjes functions are double Laplace 
transforms, this inclusion is easy to verify. For every / G § the fractional power /" G S C CM, < a < 1, is again 
a Stieltjes function. Thus, /" is completely monotone and therefore infinitely divisible. The formula 



+ 7[0,oo[ 



shows that a;~^(l + is completely monotone; but it cannot be a Stieltjes function, since it has poles at ±z, 

and ([2]) indicates that a Stieltjes function has a holomorphic extension to the cut plane C\] — cx),0]. 

Definition 4. A function / : ]0, oo[ — * [0, oo[ is called a Bernstein function, if it is of class C°° such that /' G CM. 

The set of Bernstein functions is denoted 'B3'; it is a convex cone closed under pointwise convergence. Since 
a Bernstein function is non-negative and increasing, it has a non-negative limit f{0+). Integrating the Bernstein 
representation of the completely monotone function /' gives the following integral representation of / G 23 J 

fix)=ax + P+ [ (l-e-"*)z.(di), (3) 

J]0,oo[ 

where a,/? > are constants and is the Levy measure, i.e. a positive measure on ]0, cx)[ satisfying 

/ iy{dt) < CO. 

J]0,oo[ 1 + * 

The following composition result will be useful throughout the paper, see [3j: 
Theorem 5. Let X be any of the sets TiJ, CM. Then 

f eX, g em =^ f o g eT. 
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If we assume that the representing measure iy{dt) in ^ is of the form i^{dt) = m(t) dt where m{t) is completely 
monotone, we get the family of complete Bernstein functions. We denote the collection of all complete Bernstein 
functions by CSJ. It is not hard to see that C'B?' is, Hke 23^, a convex cone which is closed under pointwise 
limits. Complete Bernstein functions are widely used in various fields and they are known under a variety of 
different names: Bondesson r2-class, operator-monotone functions, NevanHnna-Pick functions are some synonyms. 
An extensive survey can be found in SchilHng, Song and Vondracek [23], a short introduction is in [25 and [15] . 
Among the most prominent examples of complete Bernstein functions are 

Ax _ „, . , ^ 1 



X I 



(A > 0), X ^ x°' (0 < a < 1), a; i-^ log(l + x), x \/x aician ■ 



X + x 

Further examples are given below in Table [31 Many Bernstein functions given in closed form are already in C'B?'. 
The only known examples from TiJ \ QTiJ are of the form ^ Pi (1 — e"'^*'^) with pi,\i > 0. 

finite 

The following results are adapted from [23] . 

Theorem 6. A function f :]0, oo[^ [0, oo[ such that /(0+) exists is a complete Bernstein function if, and only if, 
it has an analytic extension to the cut complex plane C\] — oo,0] such that Imz ■ Im/(z) > 0, that is, f preserves 
upper and lower half planes. 

A proof can also be found in [15]. This geometric description can be used to show that, for any / ^ 0, 

/ e em 4=^ x^ ^ e § ^ x^ e em ^ i e §. 

Moreover, the following composition results are immediate. Let X be either C'B?' or §, then 

/, geX =^ foge em. 

The following stability properties are less obvious, see [23j . 
Theorem 7. Let f,g,he em and / ^ 0. Then 

i) (r (x) + e em for all a e [-1, 1] \ {0}; 

ii) {f{x")+g{x'')y^" e em for alia e[-l,l]\{0}; 

iii) fix") ■ gix^-°') e em for all a G [0, 1] . 



iv) hif{x))-g[j^^]eem. 
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All proofs, nicely cleaned up, will appear in [23]; assertion (iv) was discovered by Uchiyama [27l Lemma 2.1], 
and since fractional powers f{x) = x", < a < 1 are in CBS^, (iv) entails (iii). For positive a > assertions (i), 
(ii) are in Nakamura [20], his proofs are easily adapted to a < since / G CS?" if, and only if, 1// G §. 

Letting a ^ in Theorem [7] proves 

lirn(ir + ig")^/"^/f^ 
and since pointwise limits of complete Bernstein functions are complete Bernstein, we conclude that 

Using that ]0, 1[3 a /" is continuous, this can be used to give a new proof of the so-called log-convexity of the 
cone e'BT. 

f,ge eSJ, a e [0, 1] =^ r- g'-" G CBJ. 
An alternative proof can be found in Berg [2] and [23] . 

3 Variograms and their stability properties 

As already emphasized in Section 1, the starting point for this work is a result in Ma \T7\ Theorem 3(i)], which is 
reported below with a short alternative proof. 

Theorem 8 (Ma 2007). Let 7 : M'' ^ M be a homogeneous function. Then 

(l _ e-^Mi)^ _ e^''^'^^^)) , (4) 

Ui > 0, i = 1, 2, is a variogram if, and only if, 7(^) — \A£^\, for \ ■ \ the Euclidean norm and A a d x d matrix. 

It is natural to ask wether Ma's theorem works only for the exponential class of variograms, or wether it can 
be generalized. The subsequent result gives an answer to this problem, offering a wide class of variograms closed 
under Schur product. The result presented herein is not a mere generaHzation of Ma's result, but the key to a 
simple proof of Theorem [8] as well. 

Theorem 9. LetM.'^ 3 £, liiO ^ radially symmetric variograms for alld^N, i ~ 1,2. Let a, (3 e [0, 1] 

such that a + f3 < I and let A be a d x d matrix. Then 

f^Ao ■■=9imn92mn 

is still a variogram on M.^ for all d G N. 
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Proof. The proof of this result follows the arguments in Schilling, Song and Vondracek [23]. Indeed, if 7^ are 
radially symmetric and permissible on any d-dimensional Euclidean space, then by a theorem due to Schoenberg 
they must be of the form 7i(^) := where £ see Schoenberg [25], Bochner [Sj Theorem 4.4.4] or 

Jacob and SchiUing [TU Appendix] for alternative proofs. 

Set ha,p{x) := giix") ■ g2{x^), a; > 0. It is enough to show that h'^jj e CM. Indeed, this proves that ha,p is a 
Bernstein function and that fa,p{£.) = /iq, /3(|A^P) is a variogram. 

Elementary algebra shows that 



h{x)=x ( ) + Pgii^'^) 



Since gi S "33^, we have that g'^ e CM and x^^gi{x) G CM. So will be the compositions g'lix") and g'2{x^), gi{x°')/x" 
and g2{x^)/x^ by a straightforward application of Theorem [H Moreover, for a + P < 1, x 1-^ x"''^^^^ is completely 
monotone. The proof is finished since completely monotone functions are a convex cone closed under sum and 
product. □ 

The proof of the result above offers a considerably easier way to show Ma's result: 

Proof of Theorem[M If liO = Theorem [9] with gi{x) = 1 — cxp(— a^a;), i = 1,2, and a ^ f3 ~ ^ shows that 
(m is a variogram. 

Assume now that ([H is a variogram. Then 

e^lL^ >- -^^^ 7^(0 

is a variogram for all ai, 02 > 0, and so is the pointwise limit as ai and then 02 — ^ 0; thus, 7^(^) is a real-valued 
variogram. As such, it has a Levy-Khinchine representation 

7'(0-Q^-C+/ {l-cos{x-0)Hdx), 

Jx^O 

where Q G R'^^'^ is positive semi-definite and is a measure with |a;p Al iy{dx) < 00. Since 7(^) is homogeneous, 
we get 

7 C = — 2 ' ■ ^ = \ vQ^ . 

for the uniquely determined, positive semidefinite square root A = ^/Q of Q. □ 

Obviously, if a = /3 = 1/2 and if we take two exponential variograms we get Ma's example, see above. Thus, 
Theorem [9] enlarges considerably the scenario for the permissibility of the product of two variograms. Several 
examples of Bernstein functions may be found in [3] or in [16]; an extensive list will be included in the monograph 
by [23]. Three celebrated classes of Bernstein functions are well known in the spatial statistics literature: 



1. the Matern class ([T^) fa,^ = 1 — 2^ /r{iy){ay/xy Ki,{a\/x), x > 0, for a,u > Q and where Ky is the modified 

Bessel function of the second kind of order v\ 

2. the Cauchy class ([l3]) fa,f}{x) 1 - (1 + a:")"'', x > 0, where < a < 1 and < /3; 



Let us mention a few more stability properties that make some classes of functions appealing for its use in 
spatial statistics. We work again under the framework of radially symmetric functions. Recall that all variograms 
7 which are radially symmetric and permissible in all dimensions d = 1,2,... are of the form 7(^) = with a 

Bernstein function / G 233^. 

Proposition 10. Let j : Mf^ M. be radially symmetric for all dimensions d = 1,2, . . . such that 7(^) = for 
some g e CB^. Then 

1^12 

is a radially symmetric variogram which is permissible for every d e N. 
ii) M'^ 9 ^ i-^ , , . and £, ICPfffiTra-) are radially symmetric variograms which are permissible for every 

deN. 

Proof. Part i) can be easily proved by noting that g G G'BJ' if and only if g{x)/x is a Stieltjes transform. Since the 
reciprocal of a Stieltjes transform is completely Bernstein, we get the assertion. 

Part ii) follows straight by noting that, for g G QHJ, the function x ^ 1/ g{\/x) is a composition of the type 
a o g o a{x), where a is the Stieltjes transform a; i— > i. Since the composition cr o 17 is a Stieltjes transform, and 
since the composition of two Stieltjes transforms is in CSiF, we have the first assertion of Part ii). If we apply Part 
i) to this variogram, the second assertion follows. □ 

For further (stability) properties of the class 623?' the reader is referred to [23]; some examples of complete 
Bernstein functions are given below. 

Another interesting problem arises when quasi arithmetic operators in the sense of Hardy, Littlewood and 
Polya [H] are applied to variograms. This means that we seek conditions which preserve the permissibility of the 
underlying structure. This has been considered in [21] for quasi arithmetic composition of covariance functions. 
We believe that the same question in connection with variograms is even more challenging from the mathematical 
point of view and is equally important as far as statistics are concerned. 

Recall that a power mean is a mapping of the form {u,v) i— > 'ip^[u,v) := {u" + v°')^^", for {u,v) G and 



/ xP Y 

3. the Dagum class ([4]) fp,-y(x) := I j 




X > 0, where p, 7 G (0, 1). 



a G M. 
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Proposition 11. Let 7,; : M'' ^ M, i = 1,2, be radially symmetric variograms for all dimensions d gN. We write 
gi for the radial function such that = ffidCP)- 

i) //gi,52 € eSJ, then ^ ^ (7?(C) + 7^(0) « for all a e[-l,l]\ {0} a variogram. 

ii) Ifgi,g2eem, i/ien $ (^id^P") + 52(|CP")) is for all a e [-1,1]\{0} a variogram. 

iii) ^ 1-^ giil^f") 32(1^1^"^") is for alio <a<l a variogram. 

Proof. Assertion iii) is a simple consequence of Theorem[9l We want to mention at this point that for gi, .92 G G'BJ', 
the resulting rotationally symmetric variogram would be again of the form with h G G'BJ', see Theorem [7| 

iv) . Both i) and ii) follow immediately from[7|i) and ii), respectively. □ 

Finally, we combine two aspects treated separately until now: Given two or three isotropic variograms, we seek 
permissibility conditions for the product of ad hoc compositions. The proposition below results from a simple 
application of Theorem [7l 

Proposition 12. Let R'* 9 ^ i-^ 7i(C)j * ~ I72, 3, radially symmetric and isotropic variograms for all d E N and 
assume that 7i(^) = 5i(|^P) where gi G CS5'. Then 



are still permissible for all d E N and of the form /i(|CP) with some h e G'BJ'. 

Several examples of completely Bernstein functions for which spatial statisticians are not aware of can be found 
in the monograph [23]. Table [3] offers a small selection. 

4 Kernels and variograms of the Schoenberg-Levy type 

This section explores some results that may be obtained when working with kernels of the Schoenberg-Levy type. 
These kernels are extensively used in the literature and we refer to Ma [17] and the references given there. For 
£.i,S,2 S R'^, these are nonstationary covariance functions obtained from a nonnegative function g : [0, oo[^ [0,oo[ 
such that 5(0) = through the Hnear combination 



Celebrated examples are the fractional Brownian sheet [T], with g{^) = |^|", a e (0,2] . Ma [T^ shows that for a 
fixed ^0 G the function 




and 




5(161) + 5(161) -5(16 -61). 



Qo(6 :=5(IC + 6l)+5(k-6l)-25(l6) 
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Function 


Parameters Restriction 


Function 


Parameters Restriction 


1 


< a,/3 < 1 


ex — xf 1 H — 
V X/ 


f " 






' x + 1 


f xP y 
\i + xp) 


< 7, /9 < 1 


- - - log 1 
a X ^ 


'l + ^i) 
^ a/ 


a > 


x°' ~ x{l + x)°'-^ 
(1 + a;)" - X" 


< a < 1 


Fx sinh'^ \^^2x 
V 2 sinli(2x/2^) 






a > 


a;!-'' e°" r(z/; ax) 


a>0, 0<i^<l 


a;(l-e-2V^) 
\/x + a 


a > 


x'' e"^/^ r 




a>0, Q <v <l 



Table 1: Examples of complete Bernstein functions. r(a; x) := t"' ^ e * is the incomplete Gamma function. 

is a covariance function provided that g{\£,\) is a variogram. Indeed, we are going to show that this is a simple 
consequence of the following more general result. 

Lemma 13. Let 7(^) he a continuous variogram, and let £,,ri arguments ofW^, d G N. Then 

7^(0 :-7(e + r;)+7(e-r/)- 27(e) 
is a continuous covariance function (in £,) and a continuous variogram (in rj). Moreover, 

cj^.io ■■= 27(^7) + 27(c) - 7(e + v)- 7(e - v) 

is a continuous variogram both in ^ and in 77. 

Proof. Recall the following elementary formula for the cosine: cos(a + b) + cos(a — b) = 2 cos a cos b. Since 7(e) has 
the following Levy-Khinchine representation, 

7(C) + / {l-cosx-OHdx) 

Jx^O 

we find 

7,(0 = Q(C + »?) • (C + »?) + Q(C - ^) • (C - ^7) - 2gc • C 
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+ / (2 cos X ■ ^ — cos X ■ {^ + rj) — cos x • — 77)) ^{dx) 

Jx^O 

= 2Qij ■ rj + / (2cosx • ^ — 2cosa; • ^cosx • 77) i'(dx) 

Jx^O 

— 2Qt] • 77 + 2 / (1 — cosx • 77) cos a; • ^ v{dx). 

Jx^Q 

This shows that j^iO is symmetric and positive definite, hence a covariance function. Now consider 
= 27(77) -7,(0 

= 2(577 ■ T] + 2 (1 — cos X ■ 77) i^{dx) — 2(577 • 77 — 2 / (1 — cos x • 77) cos x ■ ^ v^dx) 
Jx^a Jx^o 



= 2 (1 — cos a; • ry) (1 — cos x ■ ^ v{dx). 

Jx^O 

Thus, (/)r;(C) is a variogram in ^ for every fixed 77 and vice versa. The proof is completed. □ 

The previous result has some obvious generalizations. Consider, for example, continuous complex-valued func- 
tions 7 : M'' ^ C satisfying the permissibility condition (U) for all ^1, . . . , ^„ £R'^. Since these conditionally negative 
definite functions also enjoy a (complex) Levy-Khinchine representation, cf. Berg and Forst [3], an argument along 
the fines of the proof of Lemma [13] shows that for every fixed ^0 G 

<^5„ (0 := 27(0 + 2Re7(eo) - - Co) - 7(C + Co) 

is permissible and has the following Levy-Khinchine representation 

ko (C) - 2 / (1 - e^y<) (1 - cos(y • Co)) '^{dy) 

Jy^O 

where v is the Levy measure of 7. Taking everywhere real parts, we are back in the variogram case of Lemma [131 

5 Eventually constant, periodic and permissible variograms 

In this section we characterize (rotationally symmetric) variograms which are constant on some annulus. Recall 
that we call a variogram eventually constant if it is constant outside some compact set K C M."^. Variograms of 
this type are, for instance, the celebrated spherical model, see Chiles and Delfiner [6J and literature quoted there, 

3,„ 1, 



7(0 2ICI+2ICI 



2 



that is a valid isotropic variogram on and lower dimensional spaces, but not in M.'^ for d > 5, see Gneiting [llj . 
Other compactly supported variograms can be obtained by ^ 1-^ Cr{0) — (^r(C) where Cr is a compactly supported 
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covariance function that vanishes outside the ball of radius r. Other models that are suitable for this purpose are, 
for instance, the Wendland class [28] 

(l-y)^, (5) 

which is permissible whenever I > [|J + 1. 

Compactly supported covariances and eventually constant variograms are particularly requested by spatial 
statisticians because they can represent effective solutions to the analysis of massive spatial dataset. The best 
linear unbiased predictor, alias kriging, depends on the inverse of either the covariance or the variogram matrix — 
and a compactly supported matrix is sparse and thus inversion is computationally less demanding. For this reason, 
the literature requires, on the one hand, new models of compactly supported correlations with desirable features, 
and, on the other hand, inspection of the general properties of compactly supported correlations or, equivalently, 
finally constant, variograms. A thorough study can be found in Wendland [29] . Of particular importance is his 
striking result in Theorem 2, Chapter 9, which entails that a candidate function that is compactly supported and 
radially symmetric cannot be positive definite for all dimensions d. This is consequence of the fact that completely 
monotone functions admit an analytic continuation to the complex half plane {Re 2; > 0}. Let us give some 
important properties that will be stated in terms of variograms. 

Proposition 14. Let j : ^ R be a variogram. If j{y) — 7(0) for some y ^ 0, then ^{^ + y) ~ j{0 for oil 
^ e M"*, i.e. 7 is periodic with period y. 

Proof. We assume 7(0) = 0, otherwise we switch to 4>{S,) := 7(^) —7(0) and note that is still a variogram. The 
following inequality is well-known for conditionally negative definite functions, see [S] proof of Proposition 7.15] 

17(0 + 7(-y) - 7(e - < (2Re7(0 - 7(0)) (2Re7(^) - 7(0)) 

(we kept the notation used in C even though 7 = Re 7 is real valued). Since 7(0) — 0, we conclude from 
7(2/) = li^y) = that 17(0 ^ 7(C ^ y)\^ 1^ 2Re7(^) -0 = which proves periodicity: 7(^ zty) — □ 

Remark 15. With very little changes the proof of Proposition [14] extends to complex valued conditionally negative 
definite functions 7 : R'' — > C and shows that 7 is periodic. 

Proposition 16. Let 7 : R — > R 6e a variogram such that j\I = for some interval I = [a, a + e] C R, a G R, 
e > 0. Then we have 7 = 0. 

Proof. Recall that square roots of variograms are subadditive, cf. 0, Prop. 7.15]. For A G (0, 1], one clearly has 
\ ■ e = {a + \e) — a and thus 

V7(Ae) < V7(a + Ae) + V7(-a) < V7(a + Ae) + V7(a) = 0. 
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We conclude that 7|(0, e] =0 and in view of Proposition [14] we get 7 = 0. □ 
Let 7 : R'^ — > R be a variogram in M'^ and write £, — {£,1, ■ ■ ■ ,£,d) = X^iLi ^e^i where = (0, . . . , 0, 1, 0, . . . , 0) is 
the jth basis vector in E''. By subadditivity we get 

d 

1=1 

Since R3 £i y-^ li^i^i) is a one-dimensional variogram, Proposition [16] shows at once 

Corollary 17. Every variogram 7 : M'' ^ M such that supp7 C R'^ is compact is identically zero. 

A closer inspection of the proof of Proposition [16] reveals that it is enough to assume that the variogram 7 
is zero on sets of the type {tej : aj < t < /3j} for some constants < < l3j and all coordinate directions 
J = 1, 2, . . . , d. 

Things are different if a variogram 7 is constant on some set but not equal to 7(0). Consider, for simplicity, 
rotationally invariant variograms, i.e. variograms 7 of the form 7(^) — and assume that /|[a,5] — 7, i.e. 7 

is constant on some annulus of width b — a and inner radius a. If we know that rj 1-^ /(I'yP) is a variogram for all 
dimensions 77 e M'*, Schoenberg's result already mentioned in Section [3] tells us that / is a Bernstein function. As 
such, / is real-analytic, and f\[a, b] can only be constant if it is constant on [0, 00). This already proves 

Proposition 18. Let 7 : M'^ R 6e a radially symmetric variogram which is permissible on any d- dimensional 
Euclidean space. Further, assume that 7(0) = 0. Then, there does not exist any ^0 > such that 7(^) — c for 
all lei >eo. 

Often we are in a situation 7(^) — is a variogram up to some dimension only. Typically such examples 

are constructed with the help of Polya's Theorem [21 as for instance the Wendland function in equation |[5|). The 
following theorem shows that this is more or less the only situation that can appear. 

Theorem 19. Let 7 : R'* ^ K, > 1, be a rotationally symmetric, continuous variogram, i.e. = for 
some f : [0, 00) [0,co). Then: 

i) R'^ 9 77 1-^ /(hp) is a variogram for all dimensions k = 1,2, . . . ,d; 

ii) / is increasing, continuous, subadditive and concave; 

iii) If f\[a,b] = c on some interval [a,b] C [0, 00), then /|[a, 00) = c. 

Proof. We write 77 for vectors in R"*^^ and ^ for vectors in R''; by (77, 0) we denote a vector in R''^^ x R = R'^. 
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Assertion i) follows easily from the observation that for rj e K*^^^ /(I^P) — /(l(^iO)P) and a straightforward 
application of Theorem [1] with £,j = (?7j,0) E M''. 

Assertion iii) follows at once from ii) since / is increasing and concave. 
To see ii) we take any 771, 77,1 S M''"^ and /ii, . . . , ^„ G C. Set for a e M 

n 

We claim that a 1-^ 0(0?) is a one-dimensional variogram. We verify the conditions of Theorem [TJ Obviously 
(f)(0) = and (/'(a) = 0(— a). Now take ai, . . . , q;„i G M and Ai, . . . , A„i S C with J2™=i '^p — 0- We find 

p,g— 1 1 

771 n 

p,q=l j,k=l 

= E E /(lfe'"p) " (Mj^p)(AifcAg) < 0. 

In the last step we used that is a variogram in as well as the characterization of variograms by Theorem 

m 

This shows that (j){a) > 0(0) for all a > 0. Thus + a^) > /(|^P) which shows that / is indeed increasing. 

That / is continuous on [0,oo) is trivial since ^ 1-^ /(ICP) is continuous. 
Moreover, 

e 5a(0 /(ICn + /(a') - /(leP + «') - /(O), e e a e M 
is a variogram. Clearly, 5a(0) = and ga{0 — ga{-~0- Choosing £,i,---,£,n G M'^ and Ai,...,A,i G C with 
= we get 

E - 6)A,A, = (/(lO - - /(le, - 6P + a')) A,A, < 

j,k=i j,k=i 

since 4>(a) > 0. Taking rt = 2, Ai=— A2=:l and = ^ and = reveals that 

/(icp + < /(o) + ./(lep + c.^) < /(len + ./(«'). 

This proves subadditivity of / on [0,oo). 
Consider the function 

0<r^h^{r) := f{r) + /(a^) _ ^ ^2)^ 

The above calculations show that ha{r) = ga{r) + /(0) is a function such that /ia(|^P) is a variogram. Thus, ha{r) 
is increasing: ha{r + c) > /iq(c) and we find with = c' and for all r > 

/(r + c) + fid) - /(r + c + c') > /(c) + /(c') ~ /(c' + c) \ (/(c) + /(2c' + c)) < /(c' + c). 
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Since c' + c = ^{c + 2c' + c) this is just the concavity of / at mid-points; as / is continuous, we get that / is 
concave. □ 

We conclude this section by another curious way to construct continuous variograms and, more generally, 
complex-valued conditionally negative definite functions with the help of Bernstein functions. The interesting fact 
in the example below is the product structure which is quite unusual for conditionally negative definite functions. 

Proposition 20. Let f be a Bernstein function such that the representing measure v in the Levy-Khinchine formula 
([3]) has a monotone decreasing density m, i.e. f{x) ^ a + (3x + oo[(-'^ ~ e~^*) m{t) dt. 

Then ^ ^ —i^f{iO is conditionally negative definite, and ^ i-^ ^ Re(*^/(*0) continuous variogram. 

Proof. By the monotonicity of m we see that m(t) = iy[t, oo) for a (Levy-) measure i.e. a measure v on ]0, oo[ 
satisfying jjp oo[*(^ ^ *)~^ iy{dt) The integration properties of ly become clear from the calculation below since we 
have only used Fubini's theorem for positive integrands to swap integrals. For a; > we get 

/"OO />oo 

xf{x)=ax + f3x^ + J a;(l-e-^*) J v{ds)dt 
= ax + f3x'^ + J x(l - e"^*) v{ds) dt 



ax 



+ (3x^+ / x{l-e-'=*)dtiy(ds) 
Jo Jo 

= ax + I3x^ + / [e"^"* - 1 sx] v{ds) 
Jo 

which, as by-product, shows that A s fi{ds) < oo. We may, therefore, plug in z = and get 

*?/(*0 = -{-ta^ + be + [1 - c-*^« - is(\ fi{ds)) . 
Thus, 7(^) := —i^f{iS,) is conditionally negative definite and Re7(^) is a variogram. □ 
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